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Abstract 



We study the eigenvalues S n £ of the Salpeter Hamiltonian H = f3\J m 2 + p 2 + vr 2 , 
v > 0, (3 > 0, in three dimensions. By using geometrical arguments we show that, 
for suitable values of P, here provided, the simple semi-classical formula 

£ n e w min \ v(P n g/r) 2 + /3a/ m 2 + r 2 \ 

provides both upper and lower energy bounds for all the eigenvalues of the problem. 
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1. Introduction 

The Bethe-Salpeter formalism [1] is generally accepted, in principle, as the appropri- 
ate framework for the description of bound states within a (relativistic) quantum 
field theory. Unfortunately, almost all applications of this formalism face serious 
problems of both conceptual and practical nature. In particular, it turns out to 
be a highly nontrivial task to extract exact information about the solutions. Con- 
sequently, one is led to consider some reasonable simplifications of the full Bethe- 
Salpeter equation, such as the following: the elimination of any dependence on time- 
like variables by the use of a static or instantaneous approximation to the interaction 
kernel, leading to the 'Salpeter equation' [2]; the neglect of the spin degrees of free- 
dom of the bound-state constituents; and the restriction to positive-energy solutions. 
With these constraints one arrives at the 'spinless Salpeter equation', which may 
be regarded as a straightforward generalization of the nonrelativistic Schrodinger 
formalism towards relativistic kinematics: it describes the bound states of scalar 
particles as well as the spin-averaged spectra of the bound states of fermions. 
In this paper we study the Salpeter Hamiltonian 

H = (3^m 2 + p 2 + vr 2 (1.1) 

in which (3 > is a parameter (allowing, for example, for more than one particle) , 
m is the mass, and vr 2 is the harmonic-oscillator potential with coupling v > 0. 
In the momentum-space representation [3-7] the operator p becomes a c-variable 
and thus, from the spectral point of view, the Hamiltonian H = j3^m 2 + p 2 + vr 2 
is equivalent to the Schrodinger operator H given by 

H = -vA + V(r), V(r) =pVm 2 + r 2 . (1.2) 

Since the potential V(r) in (1.2) increases without bound, we know [8] that the 
spectrum of the operator H is entirely discrete, and we denote its eigenvalues (also 
the eigenvalues of H) by 

£ ne (v,P,m), n = l, 2,3,..., £ = 0,1,2,..., (1.3) 

where n 'counts' the radial states in each angular-momentum subspace, labelled 
by I. Because V(r) is at once a concave function of r 2 and a convex function of 
r, this allows us to derive in Sec. 2 the approximation formula 

£ n i « min < v ( — + (3\/m 2 + r 2 > , (1.4) 
r>0 V r J 1 
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which, for suitable values of P, provides both upper and lower bounds to the energy. 
2. The energy bounds 

If we think of the potential V(r) as a smooth transformation V(r) = g(sgn(q)r q ) 
of a pure attractive power potential sgn(g)r 9 , then the potential V(r) clearly has 
two such representations, each with definite convexity. That is to say, we may write 
the potential in the following way: 



V(r) = py/m 2 + r 2 = g {1 \r 2 ) = g {2 \r), (2.1) 

where the two transformation functions g have the properties that g^ is concave 
( g" < ) and g^ = V is convex ( g" > ) . These are precisely the conditions 
under which our 'envelope theory' [9-11] applies. The tangent lines to the transfor- 
mation functions g are, on the one hand, shifted oscillator potentials of the form 
a + br 2 , and, on the other, shifted linear potentials of the form a + br. The potential 
V(r) itself is at once the envelope of an upper oscillator family and a lower linear 
family of potentials. It follows [11] that we may write 



mm < v 

r>0 V r 



Pni{l) 



) + ;\ //,-' • r- j < S nt < mm jr (^y^j + \ m- ■ /-' j , 

(2.2) 

where the 'upper' and 'lower' P numbers are defined in terms of the correspond- 
ing exact power eigenvalues E as follows. Suppose the exact eigenvalues of the 
Schrodinger operator for the pure-power potential —A + sgn((/)r 9 , q > — 1, are 
written E n i(q), then the corresponding P numbers are defined [11] by 



Pni(q) = \E ne (q)\ i2+q)/2q 



2 


1/9 


\q\ 


1/2 


_2 + q_ 




2 + q_ 





, q^O. (2.3) 



The limiting case q — > corresponds exactly to the log(r) potential, but that is 
another story [12]. Thus we have for the harmonic-oscillator potential ( q = 2 ) 

E n£ (2) = 4n + 2£ - 1, P ne {2) = 2n + £ - ^, (2.4) 
and for the linear potential ( q = 1 ) 

P ne (l) = 2(^-) 3/ \ (2.5) 
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which are listed here in Table 1. Thus we have established the principal claim of this 
paper. The minimizations in (2.2) can be carried out exactly, yielding the following 
quartic equation for r 2 : 

r 8 = ^-(m 2 + r 2 ). (2.6) 

However, given the ease of contemporary computing, the minimizations in (2.2) 
are perhaps preferable to the exact analytical solutions to (2.6); moreover, (2.2) in- 
volves meaningful semi-classical energy expressions which exhibit clearly the kinetic 
and potential energy contributions and how they depend on the parameters of the 
problem. It turns out that these energy bounds are quantitatively equivalent to the 
bounds obtained in the Appendix of Reference [7] by the use of optimized opera- 
tor inequalities. The present formulation of the energy bounds, based on convexity, 
allows for the uniform and succinct expression of our upper and lower results by 
Eq. (2.2), and admits further natural generalizations. 

As an illustration of the accuracy of (2.2) we have plotted graphs of the ap- 
proximate eigenvalues given by (2.2) in Fig. 1. We illustrate the case (3 = v = 1 and 
plot the eigenvalue bounds obtained as functions of the mass m. Corresponding 
accurate numerical eigenvalues are shown as a dashed curve between each pair of 
bound curves. These graphs confirm numerically what is immediately clear directly 
from the operator H, that the corresponding (Schrodinger) problem 7i approaches 
the oscillator for large m, and the linear potential for small m. 

3. Conclusion 

The spinless-Salpeter eigenvalue problem is not easy to solve. Even for the harmonic- 
oscillator potential, one obtains an equivalent Schrodinger problem with the poten- 
tial V(r) = f3\fm 2 + r 2 that does not admit exact analytical solutions in terms of 
known special functions. Thus, even for the oscillator problem, we must resort to 
approximations of some sort. The potential V(r) increases monotonically to infin- 
ity and we therefore know a priori that the spectrum is entirely discrete. Hence, 
the spectrum and wave functions can be found numerically with considerable ease. 
In spite of this, it is always desirable to have an eigenvalue formula, even an ap- 
proximate one, which tells us how the spectrum depends on all the parameters of 
the problem. In this paper we have used geometrical envelope theory to generate 
simple semi-classical expressions that provide upper and lower bounds to all the 
eigenvalues. 
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Table 1 Numerical values for the P numbers for the linear potential ( q = 1 ) 
used in the Schrodinger eigenvalue formula (2.3) 



n 






1 





1.37608 


2 





3.18131 


3 





4.99255 


4 





6.80514 


5 





8.61823 


1 


1 


2.37192 


2 


1 


4.15501 


3 


1 


5.95300 


4 


1 


7.75701 


5 


1 


9.56408 


1 


2 


3.37018 


2 


2 


5.14135 


3 


2 


6.92911 


4 


2 


8.72515 


5 


2 


10.52596 


1 


3 


4.36923 


2 


3 


6.13298 


3 


3 


7.91304 


4 


3 


9.70236 


5 


3 


11.49748 


1 


4 


5.36863 


2 


4 


7.12732 


3 


4 


8.90148 


4 


4 


10.68521 


5 


4 


12.47532 
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Figure 1. 

Upper and lower bounds on the energy levels of the spinless Salpeter Hamiltonian 
H with a harmonic-oscillator potential, H = f3yj m 2 + p 2 + vr 2 , for (3 = v = 1, 
as functions of the mass m. Accurate numerical eigenvalues are shown as dashed 
curves between each pair of bounds. 



